Abstract. In this paper we present calculations of pe- 
riod distribution for old accreting isolated neutron stars 
(INSs). After few billion years of evolution low velocity 
INSs come to the stage of accretion. At this stage INS's 
period evolution is governed by magnetic braking and an- 
gular momentum accreted. Since the interstellar medium 
is turbulized accreted momentum can either accelerate or 
decelerate spin of an INS, therefore the evolution of period 
has chaotic character. Our calculations show that in the 
" ' case of constant magnetic field accreting INSs have rela- 
. tively long spin periods (some hours and more, depending 
] on INS's spatial velocity, its magnetic field and density of 
■ the surrounding medium). Such periods are much longer 
bJO' than the values measured by i?OSAT for 3 radio-silent iso- 
^ . lated neutron stars. Due to long periods INSs should have 
' high spin up/down rates, p, which should fluctuate on a 
I i; time scale of ^ 1 yr. 
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1. Introduction 

Spin period is one of the most precisely determined pa- 
rameters of a neutron star (NS). Estimates of masses (for 
isolated objects), radii, temperatures, magnetic fields etc. 
are usually model dependent. That is why it is very im- 
portant to have a good picture of evolution of the only 
model independent physical parameters of INSs — spin 
period and its derivative, as they are usually used to de- 
termine other characteristics of NSs. Our aim is to obtain 
distribution of periods p for old accreting INSs (AINSs), 
we also briefly discuss possible values of period derivative 

A lot has been done to understand period evolution 
of radio pulsars (see Beskin et al. 1993) and NSs in close 
binaries ( phosh fc Lamb 1979| , [Lipunov 1992| ). AINS are 
especially interesting from the point of view of period (and 
by the way magnetic field) evolution since their history is 
not "polluted" by huge accretion, as it happens with their 
relatives in close binary systems, where NSs can accrete 
up to 1 Mq during extensive mass transfer. 

In early 90*^ there has been a great enthusiasm about 
the possibility to observe a huge population of AINSs with 
ROSAT ([Treves fc Colpi 199l|, [Blaes fc Rajagopal 1991, 



Blaes fc Madau 1993| , pVIadau 
it has become clear that AINSs 



Blaes 1994) 



However, 
very elu- 



are 

sive ([Treves et al. 199^ ) due to high spatial velocities 
or/and magnetic field properties 



( Popov et al. 2000a|) 

( Icolpi et al. 1998[ , [Livio et al. 1998[ ) 

StiU, AINSs (and radio-silent 
eral) are a subje ct of interest in 
Caraveo et al. 1996[ , [Treves et al . 2000 



INSs in gen- 
astrophysics (see 
and references 



t herein). A fe w candidates seem to be observed by ROSAT 
(Motch 2001), although it is also possible that these 



by young cooling NSs (see 


Neiihauser fc Triimper 1999 


Yakovlev et al. 1999, Walter 2001 




Popov et al. 2000bD 



Nevertheless, such sources should be very abundant 
at low fluxes available for Chandra and XMM-Newton 



observatories (in Popov et al. 2000b the authors obtain 
that at ~ 10"^'^ erg s~^ AINSs become more abundant 
than young cooling NSs, and expected number is about 
1 source per square degree for fluxes > 10~^^ erg s~^), 
so the calculation of properties of AINSs is of great 
importance now. 

Previous estimates of spin properties of AINSs 



( Lipunov fc Popov 1995a Konenkov fc Popov 1997 ) have 



given only typical values of periods, no realistic distribu- 
tions of this parameter are calculated. "Spin equilibrium" 
of NSs with the interstellar medium (ISM) has been al- 
ways assumed, i.e. the authors have considered the situ- 
ation, when all INS have enough time for spin evolution, 
they have not considered NSs with relatively high spa- 
tial velocities. In this paper we present full analysis of the 
problem. 

Period estimates are especially important as far as this 
parameter can be used to distingush accreting INSs from 
young cooling NSs and background objects. 

We proceed as follows: In the next section we describe 
the model used to calculate period distribution. In Section 
3 we present our results. In this paper we do not address 
the question of the total number of AINS, for this data 



we refer to our previous calculations (Popov et al. 2000a, 
Popov et al. 2000b). Here we only show period distribu- 



Send offprint requests to: S. Popov 



tions for AINSs. In the last section we give a brief discus- 
sion, derive typical parameters of p for AINSs and sum- 
marize the paper. 

2. Model 

In this section we describe our model of spin evolution of 
AINS in turbulent ISM. We consider constant ISM density 
and isotropic Kolmogorov turbulence^ (Kolmogorov spec- 
trum is in good correspondence with most of observation 
of interstellar, i.e. intercloud, turbulence, see for example 

^ Any other more complicated model of turbulence requires 
better knowledge of the ISM structure especially on small 
scales. This information can be obtained from radio pulsar scin - 
tillation observations (see for example Smirnova et al. 1998). 



We plan to include this data in our future calculations. 
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Falgarone & Philips 199C). Orientations of turbulent cells 
at all scales are assumed to be independent. Turbulent ve- 
locities at different scales relate to each other according to 
the Kolmogorov law: 



Vt{r2) 



1/3 



Observations show (see Ruzmaikin et al. 198^ and ref- 
erences therein) that turbulent velocity at the scale of 
i?t ~ 2 • 10^0 cm ~ 70 pc is about Vt ~ 10 km s'^. The 
scale above is close to the thickness of the gas disk of the 
Galaxy, and the velocity above — to the speed of sound 
in the ISM. It corresponds to the largest cell size possi- 
ble and the fastest movements (otherwise turbulence will 
efficiently dissipate its energy in shocks). 

As an AINS moves through the ISM it can capture 
matter inside so called Bondi (or accretion, or gravita- 
tional capture) radius, Rg = IGMjv^ . Here G is the grav- 
itational constant, M is the mass of NS, v — \/vj^g + v1 , 
VNS is the spatial velocity of a NS relative to ISM, Vg is the 
speed of sound (we assume Vs = 10 km s~^, M — IAMq 
everywhere in the paper; Vg can be dependent on the lu- 



minosity of AINS (Blaes et al. 1995), but here we neglect 
it planning to iclude this dependence in our future calcu- 
lations). 

Accretion rate M at the conditions stated above is 
equal to ( |Hoyle fc Littleton 1939| , [Bondi fc Hoyle 1944| ): 



M 



QHUlpV , 



here n is the number density of the ISM, nip is the mass 
of proton. This accretion rate corresponds to luminosity: 



L = GMM/R ■ 



lO^^n (-y/lOkms^i)"^ ergs^^ 



Based on population synthesis models 
( Popov et al. 200"0b| ) we can expect, that on average 



AINS should have luminosities about 10 erg s~ . Simple 
calculations show, that most of this energy will be emitted 
in X-rays with a typical blackbody temperature about 
0.1 keV (if due to significant magnetic field accretion pro- 
ceeds onto small polar caps then the temperature would 
be higher up to 1 keV). Note, that Bondi rate is just the 
upper limit, in reality due to heating ( phvartsman 1970a , 
Shvartsman 19711 iBlaes et al. 1995) and magnetospheric 



effects (Toropina et al. 2001) the accretion rate can be 
lower. 

Due to the turbulence accreted matter carries non-zero 
angular momentum: 



Jt 



vtiRa) ■ Rg = vtiRt)Rt ^'^R^ ■ 



In this formula it is considered that cells of the size r = Rg 
are the most important, otherwise in the Eq.(|l|) below it 
is necessary to introduce factor a ^ I. 

If it is larger than the Keplerian value on the mag- 
netosphere boundary (i.e. on the Alfven radius Ra — 



{ix"^ /2My/GMf/\ here /i is the magnetic moment of a 
NS) an accretion disk is formed around the AINS. In the 
disk a part of the angular momentum is carried outwards, 
and the NS accrets matter with Keplerian angular mo- 
mentum jK = vk{Ra) • Ra, vk - Keplerian velocity. This 
situation holds only for very low magnetic fields and low 
spatial velocities of NSs, thus we do not take it into ac- 
count in the present calculations. 

We consider the lowest spatial velocity of NSs^ to 
be equal to 10 km s~^. This value is of order of the 
speed of sound. Therefore lower spatial velocities can 
not change the accretion rate significantly. Moreover such 
low values are not very probable due to non-zero spa- 



tial velocities of NSs pro genitors (see Popov et al. 2000a, 



Arzoumanian et al. 2001 for limits onto the fraction of low 



velocity INSs derived by different methods). For several 
cases below we consider v ~ wjvs- 

During the time required to cross a turbulent cell of 
the size Rg, At = Rg/vns - 11.8 (u/lOkm s"^)-^ yrs, 
the change in the angular momentum of a NS J is given 
as: 



|AJ| = MjtAt = T:nmpVt{Rt)Rt ' R, 



l/3nl3/3 



G 



1.23 • lO^^gcm^ s" 



vtiRt 



lOkms" 



Rt 



2 • 102" cm 



-1/3 



lOkms 



-26/3 



Accordingly, the change of the spin frequency reads: 
|Aa;| = \AJ\/I 



1.23- 10"" s" 



vtjRt) 
lOkms" 



Rt 



2 • 1020 cm 



-1/3 



lOkms 



-26/3 



M5 



here / — 145 lO''^ g cm^ — moment of inertia of a NS. 
Orientation of Au; is random, and it is isotropically dis- 
tributed on the sphere. The value of the frequency change 
is strongly dependent on the spatial velocity of the NS: 
Alj ^ u~26/3^ gQ i-jj^g maximum value for ujv5 = 10 km 
is Awmax — 6 • 10~* rad s"^. 

On this view, it follows that it is possible to describe 
the spin evolution of an AINS as random movement in 3-D 
space of angular velocities u). Since typical temporal and 
"spatial" scales At and Aw are reasonably small (At <C 



to 



10^" yrs. Aw < w 



10- 



10-7 rad s" 



consider the problem to be continuous. Therefore, we can 
use differential equations valid for continuous processes. 

In this case the spin evolution of an AINS in the space 
u; is described by the diffusion equation with the coeffi- 
cient, namely: 
aAw^ 

D^^^, (1) 



6 At 



The lowest measured velocity 



^ 50 km s" 



(Lyne & Lorimer 1994), please bear in mind that selection ef- 



fects are very important here. 



4 



M.E. Prokhorov et al.: Periods of accreting isolated neutron stars 



here a — coefficient which takes into account cells with 
sizes r ^ Rq (everywhere in this paper we use a = 1). 

Besides random turbulent influence an AINS is spin- 
ning down due to magnetic braking: 

(2) 



dt "''iRljLul 
Here Rco = {GM/uj^y/^ — corotation radius, Kt — con- 
stant of the order of unity, which takes into account de- 
tails of interaction between the magnetosphere and ac- 
creted matter (everywhere in this paper we accept Ht = l)j 
Ft — random (turbulent) force with zero average value: 

On large time scale and for relatively short spin peri- 
ods we can neglect momentum of the accreted matter and 
initial period of NS po- In that case solution of Eq. (||) 
looks as: 

-1 



U) oct 



27r 



oc t . 



Magnetic braking produces a convective term in the 
evolutionary equation for spin frequency distribution, 
/{lo). As far as initial distribution of spin vectors is 
isotropic (magnetic braking does not change the orien- 
tation of u}) and turbulent diffusion is isotropic too, we 
obtain spherically symmetric distribution in the space of 
angular velocities, i.e.: 

/(u;) = Mcj) 



(/3(ijj) describes 3-D distribution of AINSs, its dimension 
From (m) and (0) wc obtain the evolutionary equation: 

here A = (fi'^/IGM). 

Boundary condition at t^; = can be derived from the 
equality of the flow of the particles at this point to zero: 



dh 
duj 



0. 



j=0 



It is easy to find stationary solution of the equation 
on the semi-infinity axi^ {u > 0) 



= Cst exp 



(4) 



3 IGMD ' 

here normalization constant Cst can be derived from the 
condition An J f^*^{uj)u!'^dLU — N, N — full number of 
AINSs in the distribution. Total number of NSs in the 
Galaxy is uncertain, N ~ 10^-10^. Population synthe- 



sis calculations (Popov et al. 2000b) give arguments for 
higher total number about lO'^. Here we do not address 
this question. From the curves below, which represent rel- 
ative period distribution of AINSs, one can determine ab- 
solute numbers of AINSs with each period value if the 
total number of these objects is known. 



^ As far as NS enter the stage of accretion with relatively 
short periods ua S> {uj{f^^)}, the obtained distribution is dif- 
fered from the real one only on high frequencies u ~ uja- 




Fig. 1. Evolution of the period distribution in time; /i = 
10'^*' G cm'^, n — 1 cm^'^, vns — 10 km s^^. Curves 
are plotted for 4 different moments from 1.72 • 10^ yrs to 
9.8-10^ yrs {tA for chosen parameters is equal to ~ 1.7-10^ 
yrs). An AINS crosses the horizontal part from ~ 10^ s to 
10"* s in ~ 6 • 10^ yrs. Curves were normalized to 1 in the 
maximux for the highest curve. 

Position of the maximum of f^^ip) depends on v, fi, 
n and M . Increasing of v and /x shifts the maximum to 
larger p, increasing of n and M — to shorter p. 

If we have to solve the problem for constant starfor- 
mation rate we can write second boundary condition as: 



IGM N 



ATTKtP- 



here N — rate at which NSs come to the 
stage of accretion, uja — '^t^/pa — accretor fre- 
quency, at that value for given v and n accre- 
tion sets on, PA = 2^/i47r(G'Af)~5/7(^2^^-)3/7 ^ 

300 yLt3/^(u/10kms-i)9/7„-3/7s. 

The function fsiw) is connected with distributions of 
absolute values of the angular velocity, fi{u!), and of the 
spin period, f{p), according to the formulae: 

32^ (27: 
fiP) = —/a — 
P \ P 

Figure 1 shows the evolution of f{p) as a function of 
time for typical parameters of an AINS and the ISM. 

For given v, p. and n NSs enter the accretor stage when 
they spin down to pA- Up to this value they evolve as 



ejectors and propellers (Lipunov 1992, Colpi et al. 2001). 

For PA < p < Per influence of the turbulence is small 
and AINS spin down according to p oc i. Here p^r defines 
the stage when spin changes due to magnetic braking and 
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turbulent acceleration/decelaration become of the same 
order of magnitude (see below the Discussion). One can 
introduce another useful timescal, Af^r- This is the time 
required to reach turbulent regime, p = Pcr- These consid- 
erations together with constant starformation rate deter- 
mine left part of the distribution where f{p) = const (see 
Fig. 1). 

As the period grows turbulence becomes more and 
more important. Finally the first AINSs reach the point 
p = oo (lu — 0) and there equilibrium component of the 
distribution (described by Eq. (^) is formed quickly. This 
component decreases in a power law fashion (/(p) oc p^^) 
at p ^ Pturb (where pturb is determined by the width of 
the distribution (^), and decreases even faster than ex- 
ponentially at p ^ Pturb- Please note that the number of 
AINSs reaching equilibrium grows linearly with time. As 
a result, the amplitude of equilibrium component of the 
distribution grows likewise. If the time required to reach 
the accretor stage, t^, is large {tgai — tA ^ tgai) the equi- 
librium component is not formed. For Ia > tgai NS never 
reach the stage of accretion. 

Similar considerations can be applied also to binary X- 



ray pulsar (see Lipunov 1992). In binaries NSs can reach 
real period equilibrium (Ghosh & Lamb 1979), and sta- 



tionary solution should be used. In that case observations 
of period fluctuations can help to derive physical param- 
eters of the systems, for example stellar wind velocity for 
wind- fed pulsars (Lipunov fc Popov 1995b). 



3. Calculations and results 

The main aim of this paper is to obtain period distribu- 
tion for AINS in turbulized ISM. We assume that NS are 
born with short (^ 1 s) spin periods. Magnetic moment 
distribution is taken in log-Gaussian form: 

1 ^ {logfi-logfio)^" 



fip) 



: exp 



(5) 



with log Ho = 30.06 and am = 0.32 (see Colpi et al. 2001 



and references therein for a recent discussion and data on 
magnetism in NSs) . Mag netic field is consi dered to be con- 
stant (see, for example, Urpin et al. 1996 for calculations 
of magnetic field decay and related discussion) . 

Velocity distribution (due to kick after supernova ex- 
plosion) is assumed to be MaxwcUian: 



f{vNs) = 



6 



exp - 



3 vjfg 



(6) 



with Vm = 200 km s 
(jy ~ 140 km s~^ (see for example Lyne & Lorimer 1994 



which corresponds to 



Cordes & Chernoff 199 



Lorimer et al. 19971 ^bordcs fc Chcrnoff 1998 



Hansen fc Phinney 1997 



Arzoumanian et al. 2001 for discussions on the 



and 
kick 



We use two values of the ISM number density, n = 1 
cm~'^ and n — 0.1 cm""^, as typical values for the Galactic 

disk g. 

All NS are divided into 30 groups in ^ interval from 
lO^S-e to lO^i-^ G-cm^ with the step 0.1 dec and 49 
groups in velocities from 10 to 500 km s^^ with the step 
10 km s~^. 

For each group we calculate tA, the time necessary to 
reach accretion: 

tA = ts + tp , 

here tE and ip — durations of the ejector and propeller 
stages correspondently. 

On the stage of ejection a NS spin down due to 
magneto-dipole losses: 



tE 



109n-i/2 ( / ^ ) f- 

vio kms-i/ v: 



yrs. 



10 kms-iy V1030 G-cm 

Spin down on the propeller stage ( ^hvartsman 1970b , 
Illarionov fc Sunyaev 1975) is not well known, but for 
constant field always tp < ts (Lipunov fc Popov 1995a) 
(see also results of numerical calculations in 
Toropin ct al. 1999), so we neglect tp in the rest of 



the paper, i.e. we assume tA — tE- 

For high spatial velocities after ejector an INS ap- 
pears not as a propeller, but as a so-called georotator 
(Lipunov 1992, see also Rutledge 2001). At this stage 



velocity of NSs). 



Ra > Rgi so surrounding matter does not feel gravita- 
tion of the NS, and its magnetosphere looks like magneto- 
spheres of the Earth, Jupiter etc. Partly because of that ef- 
fect we truncated our velocity distribution at 500 km s~^. 

For systems with tA < tgai = 10^° yrs we solve Eq. ^ 
on the time interval tgai — tA- For angular velocities we 
use grid with 200 cells from uj — to w^. Conservative 
implicit scheme is used. Derived distributions f{p,H,v) 
(for different groups) are summed with weights from Eqs. 
(|)and(|). 

Final distribution for n — 1 cm ^ and n ^ 0.1 cm 
are shown on Fig. 2. 

We note, that the final distribution significantly dif- 
fers from the one for a single set of parameters (n,/x,ti): 
compare Fig. 1 and Fig. 2. In the final one we see power- 
law cutoff (/ cx p^'^) at long periods {p > 10^ s) similar 
to cutoff on the Fig. 1. Then at the intermediate periods 
(lO'' < p < 10'' s) the curve has power-law ascent due to 
summation of individual curves maxima. Finally a sharp 
cutoff at short periods {p < 10^ s) is presented because dif- 
ferent INSs enter the accretor stage with different periods: 

PA = PA{n,lJL,v). 

4. Discussion and conclusions 

After an INS comes to the stage of accretion [p > pa) its 
spin is controlled by two processes (see Eq.2): magnetic 

^ This assumption is reasonable for relatively low velocity 
INS, but only this fraction of the whole population is important 
for us here. 
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Fig. 2. Period distribution for populations of AINSs evolv- 
ing in the ISM with number density n — 1 cm~^ (upper 
solid curve) and n = 0.1 cni^"^ (the second dotted curve). 
Lower curves show results for low velocity part of AINS 
population for n = 0.1 cm~'^ (v < 60, 30, 15 km corre- 
spondently, see Discussion below). Curves for normalized 
to 1 in maximum for the highest curve. 



spin-down and turbulent spin-up/spin-down. Schematic 
view of such evolution process is shown on Fig. 3. 

We can describe them with characteristic timescales: 
tmag and tturb- As One can calculate, tmag ~ P and Uurb ~ 
p~^. Here we evaluate t as p/p, and pturb = p^Mj /{2t:I). 
Here we again note, that jt can be larger that the Kep- 
lerian value of angular momentum on the magnetosphere, 
JK = vk{Ra), Ra- So we write j instead of jt as we did 
before, now j — Tii\a.{jK,jt)- 

Pmag = 2^^i''/{IGM) (7) 

Initially (immediately after p = pa) magnetic spin- 
down is more significant: 
IGM 

^P, (8) 

'cr, these two timescales should be- 



t 



but at some period, Pc 
come of the same order, and for longer periods an INS will 
be governed mainly by turbulent forces. 
One can obtain the following formula: 



Pc 



GMMj 

This period lies between pA and pt^ 



(9) 

An INS reach 



Per in Atcr ~ 10^-10^ yrs after onset of accretion: 



iVgm 



(10) 



Here Atcr are calculated approximately as Pcr/Pmag- 

The period evolution of an AINS can be described in 
the following way (Fig. 3): after the INS comes to the stage 
of accretion it spins down for ~ 10^ — 10 yrs up to Pcr, 




dt= Rg/v 



Fig. 3. Schematic view of period evolution of an INS. 
A NS starts with p — pq, then spins down according 
to magneto-dipole formula up to pE, then a short pro- 
peller stage (marked by a circle) appears and lasts down to 
p = PA- At the situation illustrated in the figure Ie ~ Ia- 
On the stage of accretion at first magnetic braking is more 
important, after at t = tcr turbulence starts to influence 
significantly spin evolution of a NS period stars to fluctu- 
ate with a typical time scale Rg/vns close to mean value 
noted Pturb on the vertical axis. 



then the evolution is mainly controlled by the turbulence, 
and period fluctuates with typical value pturb, which is 
determined by the properties of the surrounding ISM and 
spatial velocity of an INS. 

We note, that we do not take into account any selec- 
tion effects. For example, as far as period and luminosity 
both depend on the velocity of an INS, they are corre- 
lated. Taking into account flux limits of the present day 
satellites it is possible to calculate probability to observe 
an accreting INS with some period. Also periods (and lu- 
minosities) are correlated with position of an INS in the 
Galaxy. 

So, it is reasonable to make calculations which include 
all that effects in order to make better predictions for ob- 
servations. We plan to unite our population synthesis cal- 
culations with detailed calculations of spin evolution later. 
In this paper we present distributions as if all accreting 
INSs can be observed. 

For the main part of its life spin period of a low- 
velocity AINS is governed by turbulent forces. Charac- 
teristic timescale for them can be written as: tturb = 
Iuj/{Mj). And it is equal to lO'^-lO^ yrs for typical pa- 
rameters. 
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For field decay picture should be completely different 



(see for exanrple Konenkov & Popov 1997, Wang 1997), 
and observations of AINSs can put important lim- 
its onto models of magnetic field decay in NS 



(Popov & Prokhorov 200C). For decaying fields AINSs can 
appear as pulsating sources with periods about 10 s and 
p about 10"" s/s dPopov fc Konenkov 19"98[). The value 



and sign of p will fluctuate as an INS passes through the 
turbulent cell on a time scale Rq/vns, which is about a 
year for typical parameters. Irregular fluctuations of p on 
that time scale can be significant indications for accretion 
(vs. cooling) nature of observed luminosity. 

Roughly p can be estimated from the expression: 



\p\ < P 



2^1 

27r/ 



-17/3 



(11) 



Here we neglect magnetic braking in Eq. (2). This equa- 
tion for p is valid for turbulent regime of spin evolution, 
i.e. for p 3> pturb- For the given p value of p fluctuates be- 



tween +P ^^'^ ~P 27r/' ^^"-^ depends on n, but not 
on /J, (if j = it, not j = ja)- In this picture p distribution 
in the interval specified above (Eq.ll) is flat. 

Behavior of p and p of AIN Ss in molecular clouds can 
be different ( Colpi et al. 199^^ ), especially for low spatial 
velocities of NSs. Some of our assumptions in that case 
are not valid, and results cannot be applied directly. But 
we note, that passages through molecular clouds are rel- 
atively rare and short, they cannot signiflcantly influnce 
the general picture of AINSs spin evolution. 

Period distribution which can be obtained from obser- 
vation (for example from ROSAT data) can be different 
from two upper curves on the Fig. 2. Such surveys are flux- 
limited, so they include (in the case of AINS) the most 
luminous objects. But they form only a small fraction of 
the whole population. 

To illustrate it on the Fig. 2 we also plot distribu- 
tions for low velocity objects {v < 60,30,15 km s~^). 
In the case of fixed ISM density an upper limit on the 
value of the spatial velocity corresponds to a lower limit 
on the accreting luminosity of an AINS. Clearly, brighter 
the source — shorter (on average) its spin period. Even 
in that groups with relatively short periods their values 
are far from typical periods of ROSAT INSs, ~ 5 — 20 s. 
So, these objects cannot be explained by accretors with 
constant field B ~ 10^^ G. 

Calculations of period distributions for decaying mag- 
netic field, for populations of INSs with significant frac- 
tion of magnetars and for accretion rate different from the 
standart Bondi-Hiyle-Littleton value (due to heating and 
influence of a magnetosphere) will be done in a separate 
paper. 

In conclusion we stress readers attention on the main 
results of the paper: 

— we obtained spin period distributions for AINSs with 
constant magnetic field and "pulsar" properties (field, 
initial periods and velocity distributions). 



— these distributions are shown in Fig. 2. They have 
broad maximum at very long periods, ~ 10^ s. In that 
case observed objects should not show any periodicity. 

— periods of these objects should fluctuate on a time scale 
Rg/vns ~ 1 yr. 
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